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ABSTRACT 

We report the detection of two new planets from the Anglo-Australian Planet 
Search. These planets orbit two stars each previously known to host one planet. 
The new planet orbiting HD 142 has a period of 6005±427 days, and a minimum 
mass of 5.3Mj up . HD 142c is thus a new Jupiter analog: a gas-giant planet with 
a long period and low eccentricity (e = 0.21 ± 0.07). The second planet in the 
HD 159868 system has a period of 352.3±1.3 days, and m sin i=0.73±0.05 M Jup . 
In both of these systems, including the additional planets in the fitting process 
significantly reduced the eccentricity of the original planet. These systems are 
thus examples of how multiple-planet systems can masquerade as moderately 
eccentric single-planet systems. 
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Subject headings: planetary systems 
vidual (HD 142, HD 159868) 



techniques: radial velocities - stars: indi- 



1. Introduction 

Recent discoveries of multi-planet systems from maturing high-precision Doppler planet 
searches are revealing a surprising diversity of planetary system properties. It is becoming 
apparent that, given a sufficient number of high-precision observations, many seemingly 
solitary stars or single-planet systems are found to host additional orbiting bodies. This 
trend is evident both at extremely low masses (e.g. Anglada-Escude et al. 2012, Mayor et 
al. 2011, Vogt et al. 2010) and long periods (e.g. HD 134987 c, Jones et al. 2010; 47 UMa d, 
Gregory & Fischer 2010). 

Confirmation by independent observatories is extremely useful when testing the poten- 
tial detection of additional planets, particularly if they have sma ll radial-velocity a mplitudes 
or the host star has a high level of velocity jitter. For example, iBean et al.l ( 120081 ) reported 
a third planet in the HD 74156 system using data from the Hobby-Eberly Telescope. How- 
ever, using the same sp e ctra p rocessed with two different, independent Doppler velocity 
codes, IWittenmver et al. ( 20091) could not confirm that planet. The Keck observations pre- 
sented by iMeschiari et al.l ( 1201 ll ) were also inconsistent with a third planet in that system. 



High-precision and high-cadence data from mul tiple sites have proved critical in the con- 



firmation of low-mass plane ts such as HP 43 08 (lUdry et al.ll2006l; lO'Toole et al 



61 Vir three-planet system (jVogt et al.ll2010l ). and HD 114613 ( ITuomi et al.ll2012[ ). For can- 



20091 ). the 



didate multiple-planet systems, dynamical stability modelling is also a critical tool, since 
periodic signals arising from observational sampling or stellar activity can be misinterpreted 
as planets. The inclusion of rigorous dynamical modelling has recently shown some candi- 
date planetary systems to be unfeasible, e.g. Horner et al. (2011), Tuomi (2011), Wittenmyer 
et al. (2012). 

In this work, we present new data from the Anglo-Australian Planet Search (AAPS) 
which provide evidence for one additional planet orbiting HD 142 (§3.1) with a period of 
6005 days. We also perform a Bayesian analysis (§3.2), yielding results that are in agreement 
with the conventional least squares solution. In Section 3.3, we present AAPS and Keck data 
indicating a second planet in the HD 159868 system with an orbital period of 352 days. These 
two proposed multi-planet systems are subjected to detailed dynamical testing in Section 4. 
In Section 5 we use direct imaging of the HD 142 system to rule out the known stellar 
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companion as the source of the observed radial-velocity signal. Finally in Section 6, we place 
these discoveries in the context of the overall distribution of exoplanet properties and give 
our conclusions. 



Observations and Stellar Parameters 



AAPS Doppler measurements are made with the UCLES echelle spectrograph ( Diego et al 



1991 ). Keck Doppler measurements are made with the HIRES spectrograph (jVogt et al. 



19941 ) . An iodine absorption cell provides wavelength calibration from 5000 to 6200 A. The 
spectrograph point-spread function and wavelength calibration is deri ved from the i odine 



absorption lines e mbedded on every pixel of the spectrum by the cell (IValenti et al.l Il995 



Butler et al.lll9961 ). The result is a precision Doppler velocity estimate for each epoch, along 
with an internal uncertainty estimate, which includes the effects of photon-counting uncer- 
tainties, residual errors in the spectrograph PSF model, and variation in the underlying 
spectrum between the iodine-free template, and epoch spectra observed through the iodine 
cell. All velocities are measured relative to the zero-point defined by the template observa- 
tion. 

HD 142 has been observed by the AAT at 82 epochs, with a total data span of 5067 
days. HD 159868 has been observed by the AAT at 47 epochs, with a total data span of 3396 
days. In this analysis, we also add 34 Keck epochs spanning 1593 days. The radial-velocity 
data for HD 142 are presented in Table [lj and the d ata for HD 159868 a re in Tables [2] 
and |3j The physical par ameters of HP 142 ar e given in iTinney et al.l (120021 ) , and those for 
HD 159868 are given in lO'Toole et al.l (120071 ). Tables H] and [5] summarise the parameters 
for HD 142 and HD 159868, respectively. Briefly, HD 142 i s a GUV star with a mass of 
1.15±0.10 M , is possibly slightly evolved (jTinney et al.ll2002l ). is chromospherically inactive 
(logR' HK =-4.95) , with a moderately rapid rotation rate (V sin i=10A kms -1 ). HD 159868 
is a G5 dwarf with a mass of l.O87l°;o33M , solar metallicity, and is a chromospherically 
inactive (lo g.R / R ^=-4.96) slow rotator (v sin i—2.1 kms -1 ). HD 142 has a known stellar 



companion (IPoveda et a 



(lEggenberger et al. 



2007 



199 4j), which is a lat e K/early M star with a mass of ~0.56M Q 



Raghavan et al.l 120061 ) and a projected separation of 105.1 AU. 
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3. Orbit Fitting and Planetary Parameters 



3.1. HD 142 



The Jupiter-mass planet orbiti ng HD 142 with a p eriod of 339 days was one of the first 
discoveries reported by the AAPS (ITinney et al.ll2002l ). A further ten years of observations 
have revealed evidence for a long-period signal consistent with a second planet in the system. 
In the discovery paper, the one-planet fit had a root-mean-square (RMS) residual scatter of 
5. 9 ms" 1 , c onsist ent with the 3-4 ms" 1 estimated jitter for that star. Using the formulation 
of IWrightl (120051 ). we now estimate a jitter of 4.5 ms -1 , which we apply in quadrature to 
the internal uncertainties shown in Table [TJ A one-planet fit to the current data set for 
HD 142 now has an RMS of 31.1 ms -1 . This significantly exceeds both the scatter due to 
the underlying precision of our AAPS Doppler measurement system and the predicted levels 
of stellar activity jitter expected for HD 142. This excess scatter led us to investigate the 
possibility of one or more additional planets orbiting this star. 

Visual inspection of the residuals to the one-planet fit shows an obvious long-period 
signal (P > 1000 days), and the periodogram shows a significant peak at very long periods 
(Figured]), so we proceed to fit a second planet. First, we explored the vast and uncertain 
parameter space of the long-period signal with a genetic algorithm (e.g. Cochran et al. 2007, 
Tinney et al. 2011, Wittenmyer et al. 2012). We allowed the second planet to take on periods 
between 1000 and 10000 days, and an eccentricity e < 0.6. The genetic algorithm ran for 
50,000 iterations, each of which consisted of typically 1000-3000 generations, during which 
the two-planet fits evolved toward a x 2 minimum. The best-fit system parameters from 
this process are thus the result of ~ 10 s trial Keplerian fits. Used in this way, the genetic 
algorithm is an effective way of exploring a large parameter space, which is particularly 
important when the candidate planet's period is comparable to the length of th e available 
data. We then used the GaussFit least-squares fitting code (jJefferys et al.lll987l ) to obtain 
a Keplerian model fit, with the best 2-planet fit parameters from the genetic selection as 
initial inputs. 

The two-planet fit has an RMS of 11.2 ms -1 , which is still somewhat higher than that 
expected for this star. A periodogram of the residuals shows a pe ak at 108 days (Figured]). 
We used a bootstrap randomisation process (IKurster et al.lll997l ) to assess the false-alarm 
probability of the peak at 108 days. The bootstrap method randomly shuffles the velocity 
observations while keeping the times of observation fixed. The periodogram of this shuffled 
data set is then computed and its highest peak recorded. From 10,000 such realisations, the 
peak at 108 days has a bootstrap false-alarm probability of 5.1%. In mid-2011, this false- 
alarm probability was 2.5% - that the addition of new data did not improve the statistical 
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credibility of the 108-day signal leads us to conclude that the signal cannot be claimed as 
planetary in origin at this time. 

The residual scatter about the two-planet fit remains much higher than expected given 
the jitter estimate of 4.45 ms _1 for HD 142. However, the estimation of stellar activity jitter 
is a rather imprecise process, with uncertainties up to a factor of 2 (J. Wright, priv. comm.). 
Hence, it is possible that we have underestimated the activity jitter for HD 142. A jitter 
estimate of 11.3 ms" 1 is required to produce a reduced x 2 of unity for the two-planet fit. 
Examining the distribution of jitter estimates for similar stars (Wright 2005, top panel of 
their Figure 7), the distribution h as a tail extending t oward a maximum jitter of 11ms -1 , 
but only includes 36 stars in total. Ilsaacson fc Fischer! (120101 ) provide a different formulation 
to estimate acti vity jitter, using the Ca II S-in dex. The Mount Wi l son S -index for HD 142 
is Smw — 0.187 (I Jenkins et al.ll2006l ). Using the Ilsaacson fc Fischer! (120101 ) formulation, this 
yields an estimated jitter of 3.07m s -1 . Based on the chromospheric activity index and the 
high log g (Table HJ), it appears that HD 142 is an inactive dwarf star, which nonetheless 
presents a high level of radial-velocity noise. We note that this star is quite a rapid rotator: 
it has a v sini of 10.4 kms -1 , compared to typical planet-search targets which have v sini 
of 2-4 kms -1 . We conclude that the poor velocity precision for HD 142 is attributable to the 
rapid rotation, which broadens the spectral lines and limits our ability to derive extremely 
precise radial velocities. 



Using a stellar mass of 1.15±0.10 M (jTinney et al.l 120021 ). we estimate the minimum 
mass, m sin i, for the outermost planet (planet c) to be 5.3±0.7 M Jup . The 2-planet fit is 
shown in Figure [2] and the planetary parameters are given in Table [6j The individual fits for 
each of the two planets are shown in Figure [31 



3.2. A Bayesian Analysis for HD 142 

Thirteen years of AAT data have provided evidence for a very long-period planet 
(HD 142c) with a period of more than 6000 days. The baseline of these AAT observations 
is the longest currently available at high precision for this star, which makes independent 
confirmation of HD 142c problematic. A third candidate signal is present, with a period of 
108 days and a velocity semiamplitude K ~ 12ms -1 , which is comparable to the 11.2ms -1 
residual velocity scatter about the 2-planet fit. It is prudent, then, to employ an independent 
analysis to test the plausibility of the 108-day signal. 



We analysed the AAPS radial velocities using posterior samplings of different models and 
the comparisons of these models using Bayesian model probabilities (e.g. lTuomi &: Kotiranta 
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2009 



Tuoml2011 : iTuomi et al 11201 lh . We used the adaptive Metropolis algorithm ( jHaario et al.l 



20011 ) for posterior samplings because it appears t o be a reasonably efficient meth od for 
analysing radial velocities with Keplerian models (ITuomil 1201 ll ; ITuomi et al.l l201ll ). We 
present the results using the maximum a posteriori (MAP) estimates and th e corresponding 
99% c redibility intervals (D0.99), i-e. the Bayesian credibility sets as defined in lTuomi fc Kotiranta 
(2009). Our prior probability densities of the model parameters are those used in lFord fc Gregory 
(120071 ). with slight modifications. We penalised very high eccentricities by setting the prior 
densities for orbital eccentricities 7r(e) oc A/"(0, o~1), where the parameter a e was set to have 
a value of 0.3 that still allows the orbital eccentricities to have high values if the data insists 
so. We also adopted conservative prior probabilities for models with k Keplerian signals 
such that P(M. k ) = 2P(Mk+i), i-e. that the prior probability of having k + 1 planets in 
the system is always two times less than having k planets. Essentially, this enables us to 
be more confident with the interpretation of our model probabilities - if there appear to be 
k Keplerian signals in the data, we actually underestimate the significance of the weakest 
signal because of these priors. 

The probabilities of models Aik, k = 0, ...,3 planets and the corresponding RMS values 
are shown in Tabled The velocity jitter was also included as a free parameter in the Bayesian 
model, and the best-fit jitter values are given in Table[7]for each of the fc-planet models. These 
results support the presence of three Keplerian signals in the data, although the parameters 
of the 108-day signal ("planet d") are poorly constrained. The orbital parameters of this 
three-planet solution are shown in Table El The distribution of allowed orbital periods for the 
outer planet had a substantially longer tail toward longer periods (as expected given the data 
span is shorter than the expected orbital period for this planet). The large uncertainty in 
the period of planet c maps directly into large uncertainties for the eccentricity and velocity 
semiamplitude K. We therefore note that the nominal la uncertainties in the parameters of 
planet c as given in Table [6] are likely underestimated, and we advise the reader to consider 
the 99% confidence intervals in Table [8] as more comprehensive. 



3.3. HD 159868 



O' Toole et al.l (120071 ) reported the detection of a long-period (P = 986 days), eccentric 
(e = 0.69) planet orbiting HD 159868, based on 4.5 years of AAT data. That fit had a 
residual RMS scatter of 8.4 ms" 1 , which the authors noted as larger than expected for that 
star. They speculated that a second planet with P = 180 days would substantially reduce 
the RMS but, wary of the sampling difficulties in constraining planet candidates with periods 
near one-half of a sidereal year, they presented only the single-planet solution for HD 159868. 
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Now, with nearly twice as much data (Table [2]), the single-planet fit has an RMS of 15.8 
ms _1 . As is the case for HD 142 above, the worsening single-planet fit gives a clue that 
additional planets are present in this system. Here we adopt a jitter estimate of 2.65 ms 4 , 
and apply this in quadrature to the uncertainties given in Table |5] before performing orbital 
fitting. The fitting procedures used followed those outlined above for HD 142. 

Interestingly, a 1-planet fit for HD 159868 now has an eccent ricity of only 0.16=1=0 .11, 
which is markedly different from the e = 0.69 solution presented in lO' Toole et al.l ( 120071 ). A 
periodogram of the residuals to this fit (Figure 0J shows a large peak at 355 days. This peak 
has a bootstrap false-alarm probability <0.01%. Now armed with substantial evidence for 
a second planet in the system, we proceed with a 2-Keplerian solution. For the final orbit 
fit, we also include 4.3 years of data from the Keck telescope. The best fit has a second 
planet with a perio d of 352.3±1.3 days (Table E]); adopting a stellar mass of 1.087+q;q33 
( iTakeda et al.l 120071 ) . the planet has m sin i = 0.73±0.05 Mj up . Figure [5] shows the 2-planet 
fit and the phase coverage for the new planet candidate, which has a period near one year. 
This fit has a total residual RMS of 5.8 ms" 1 (AAT - 6.7 ms" 1 ; Keck - 4.6 ms _1 ), which 
is still somewhat higher than expected based on the instrumental noise and stellar activity 
jitter. Based on the somewhat low log g values given in Table [5J HD 159868 may be slightly 
evolved. If HD 159868 is a subgiant, the velocity jitter may be closer to ^ms" 1 , ty pical 



of subgiants ( IKjeldsen &: Bedding! Il995t IJohnson et al.ll2010l ; IWittenmyer et al.l 1201 lbh and 
consistent with the RMS scatter about our two-planet fit. 

A periodogram of the residuals now has a peak at 12.6 days, but it is not significant, 
with a bootstrap false-alarm probability of 27%. As a further check, we also performed a 
Bayesian analysis as described above for the HD 142 system, and we obtain model probabil- 
ities (Table |9]) which confidently indicate two Keplerian signals. These results demonstrate 
how two planets in nearly-circular orbits can mimic a single eccentric planet when data are 
sparse a nd more subject to vagaries of sampling. This c an cause signals near one year to be 
missed (jTinney et al.l 1201 ll ; lAnglada-Escude et al.ll2010l ). 



4. Dynamical Stability Analysis 



While single planet systems can be fully solved with a simple Keplerian analysis, grav- 
itationally interacting systems of multiple planets require a full Newtonian analysis. Pro- 
posed solutions need to be shown to be dynamically stable over reasonably long timescales 
( IHorner et al.l 1201 ll ; IWittenmyer et al.l |2012| ; iHorner et al.l 120121 ) to be considere d "real." 



Gravitationally interacting systems provide both an independent check on "reality" ( IFabrycky et al. 



20121 ) . and a useful means to constrain or solve for the orbital inclination angle and true mass 



of the planets ( IRivera et al.ll2005f ). 



4.1. The HD 142 System 

To examine the dynamical stability of the proposed HD 142 planetary system, we per- 
formed a series of highly detailed n-body dynamical simulations of a wide range of potential 
system architectures. Given the extreme uncertainty in the orbit of the outermost planet 
detected in the HD 142 system, we concentrated solely on the orbital stability of the 350- 
day planet and a possible 108-day planet. As discussed in Section 3.1, there is a residual 
signal after fitting two planets; the false-alarm probability of that signal is not presently low 
enough to justify claiming a third planet. However, if further observational data support 
the existence of such an object, it would be wise to understand the dynamics of the system. 
In this subsection, we explore the dynamical interactions between the known 350-day planet 
and the possible 108-day planet. 



Following previous studies of exoplanetary stability ( 


Marsha 


et al. 


2010: Horner et al. 


2011: Wittenmver et al. 2012: Horner et al. 2012: Robertson et al. 


20121). we used the Hvbrid 


integrator within the n-body dynamics package MERCURY ( 


Chambers 


1 1999) to examine 



the stability of the planetary system as a function of its orbital architecture. Following those 
earlier works, we considered two-planet systems in which planet b (the most well-constrained) 
was placed on its nominal best-fit orbit (Table [6]). The best fit for a potential third planet 
has P = 108.2±0.2 days, K = 12.3±1.9 ms" 1 , e = 0.28±0.16, and u = 271±27 degrees. 
The initial orbital elements were then uniformly distributed across the 3<r confidence range 
in a, e, and mean anomaly M. In total, we tested 35 values of a, 35 values of e, and 25 
values of M, creating a grid of 30,625 initial system architectures. The semi-major axes for 
the innermost planet varied between 0.44 and 0.50, the eccentricity was varied between 0.00 
and 0.76, and the mean anomalies were varied between 5 and 355 degrees. Each of these 
systems was then integrated for a period of 100 million years. A body was considered to be 
ejected from the system if it reached a distance of 10 AU from the central star. For each 
of the 30,625 potential systems, we obtained either the time at which the system fell apart 
(through collisions or ejections), or alternatively found that the system remained intact until 
the end of our simulations. 

These results allow us to construct detailed dynamical maps of the system. Figure [6] 
shows the results: detailed dynamical maps of the HD 142 system. Each panel shows the 
mean lifetime of the system, as a function of semi-major axis and eccentricity, with each 
coloured square showing the mean of the twenty-five individual runs carried out at that 
particular a-e location. Panel (a) shows the results using the nominal best-fit orbit for 
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planet b: the entire region spanned by the ± la uncertainties on the orbit of planet d 
is dynamically stable. The only departures from stable solutions are found at relatively 
large orbital eccentricities for planet d, which cause the two planets to experience mutually 
destabilising encounters. 

Following these test integrations, we examined the most extreme case possible within the 
3a error bounds on the orbit of planet b. We repeated the integrations exactly as described 
above, but placed planet b on the most eccentric orbit allowed within the 3a confidence 
interval, and at the smallest semi-major axis that interval would allow (i.e. a = 0.93 AU, 
e = 0.37). This was designed to give the system the greatest possible chance of instability 
- essentially to test it to destruction. These results are shown in panel (d) of Figure |6j In 
contrast to the integrations described above, the great majority of the phase space tested for 
this extreme scenario turns out to be unstable, although there remain several broad "islands 
of stability." These islands are separated by a wealth of unstable regions, primarily driven 
by the web of mutual mean-motion resonances between the two planets. 

To better illustrate how the region of instability varies as a function of the orbit of a 
potential 108-day planet "d," we carried out two subsidiary suites of integrations. These 
again covered the full ± 3a range of a-e space for planet d, but with a resolution of 25x25x9 
steps in a-e-M. First, we placed planet b on an orbit with eccentricity la greater than the 
nominal value, and a semi-major axis la smaller than the nominal value (i.e. a = 0. 99 AU, 
e = 0.27). In the second suite, planet b was placed on an orbit 2a more eccentric than, 
and 2a inside, the nominal values (i.e. a = 0. 96 AU, e = 0.32). These results are shown 
in panels (b) and (c) of Figure El respectively. Thus, panels (b)-(d) show the results when 
the orbital eccentricity and semi-major axis of planet b are changed in la steps from their 
nominal values. 

It is clear that the scenarios featured in panels (a) and (b) of Figure El reveal a far 
greater proportion of dynamically stable orbits. One obvious result from these dynamical 
tests is that when the eccentricity of a candidate planet d increases, the stability of the 
system dramatically decreases. The cross at the center of each panel in Figure [6] shows 
the best fit and la uncertainties for a third planet. The most recent data for HD 142 now 
result in a higher eccentricity for a third planet (e = 0.4±0.1). The results of our dynamical 
simulations, in which higher eccentricities for a proposed planet d were less likely to remain 
stable, are in agreement with the analysis in Section 3.1: a 108-day planet in this system is 
increasingly unlikely. 
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4.2. The HD 159868 System 

To test the orbital stability of the two planets discovered in the HD 159868 system, 
we once again performed two detailed suites of dynamical integrations of the planetary 
system, using the Hybrid integrator within MERCURY. In the first suite, as for the first set 
of runs performed for the planets in the HD 142 system, we considered scenarios in which 
HD 159868 b was placed on its nominal orbit (a = 2.25 AU, e = 0.05, etc.). We then carried 
out 30,625 individual simulations, through which the orbit of HD 159868 c was varied across 
the full 3-0" range of allowed orbital solutions in a, e and M (in a 35x35x25 grid, again as 
before). In stark contrast to the results for HD142, which featured a significant number of 
unstable solutions, every single system tested for HD 159868 remained dynamically stable 
for the full 100 Myr of our study. 

Following that first suite of integrations, we carried out a second test, in which the 
orbit of HD 159868 b was set to the most extreme allowed within the 3-cr uncertainties, with 
a = 2.13 AU and e = 0.17. Once again, we tested 30,625 unique planetary systems, with 
the orbit of HD 159868 c varied across the full 3-cr range of allowed orbital solutions. These 
extreme runs did yield a small number of unstable solutions (162 of the 30,625 runs were 
destabilised by the end of the integrations), but the vast majority of systems tested survived 
unscathed until the end of the simulations. Every unstable solution required the initial 
orbit of HD 159868 c to have an eccentricity of at least 0.32, and almost all featured initial 
semi- major axes of greater than 1.032 AU. 



5. Direct Imaging for HD 142 

Since the outermost object in the HD 142 system has a very long and poorly-constrained 
orbital period, it is prudent to check whether the observed radial-velocity variation is due 
to a stellar companion on a much longer-period orbit. As noted in Section 2, HD 142 is 
known to host a stellar companion (0.56M Q ) with a projected separation of 105.1 AU. If 
this companion were the source of the radial-velocity signature attributed to the outermost 
companion, it would require the system to be almost face on (an inclination angle of ~ 0.11 
degrees for the derived e = 0.20 and K = 55.5ms _1 ). 

Co mbining previous VLT-NA CO observations of the HD 142 system, previously pub- 



lished in lEggenberger et al.l (120071 ) , and a more recent observation we made with the Near- 
Infrared Coronographic Imager (NIC!) on the 8m Gemini Telescope (Figure [7] and Table [TU]) . 
we can clearly see the stellar companion moving almost directly towards the star which sug- 
gests a nearly edge-on system. The stellar companion seen from direct imaging would also 
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have an orbital period of > 1000 yr instead of the 6005 days (17 yr) determined from the 
radial- velocity data. This evidence leads us to believe that the stellar companion is not likely 
to be the cause of the radial-velocity variation attributed to HD 142c. 



Discussion 



In IWittenmyer et al.l (j2011al ) , we defined a Jupiter analog as a gas-giant planet with 
a period P>8 yr and a small eccentricity (e<0.2). HD 142c has a period of 17 years, 
e = 0.2, and a mass estimate consistent with a gas-giant planet. This planet thus represents 
a new Jupiter analog, the fourth such planet discovered by the A APS. The three pre vious 
AAPS Jupite r analogs are HP 1349 87c J Jones et alJboioh . G J 832b faailev et alil2009h . and 
HD 160691c ([McCarthy et al.l 12004 ). This discovery of a new, very long-period planet pro- 
vides additional evidence that continued support for the AAPS is bearing fruit. We note 
that although the best-fit period for HD 142c is lo nger than the dura t ion of observations, 
the available data cover ~80% of an orbital cycle. IWittenmyer et al.l (J2011aJ) investigated 
the extant literature and found that the minimum orbital coverage for published planets was 
70% of an orbital cycle. 

The AAPS data also show hints of a third signal with a period of 108 days, evident 
in both the traditional periodogram analysis and the Bayesian analysis (§3.2). Since this 
signal still has a false-alarm probability of 5.1%, we do not claim it to be a planet at this 
time. We note that while very low-amplitude planets have been detected, with amplitudes 
comparable to the stellar radial- velocity noise (e.g. Vogt et al. 2010, Pepe et al. 2011), those 
detections were clearly evident in the periodograms, and the host stars were slow rotators 
with extremely low intrinsic jitter (unlike HD 142). 

These discoveries also highlight the importance of continuing to monitor known plan- 
etary systems for signs of additional objects. In particular, both HD 142 and HD 159868 
had moderately eccentric orbital solutions as well as excess scatt er about the single - plane t 
fit. A dedicated search pursuing this strategy was performed by IWittenmyer et al.l (120091 ). 
who observed 22 known planetary systems for 3 years using the Hobby-Eberly Telescope. 
While that survey did not result in new plan et dscoveries, the new data cast do ubts on the 
existence of the proposed p lanets HD 20367b (jUdry et al.ll2003l ) and HD 74156d (IBean et al. 
20081 : iMeschiari et alJlioilh . 



In both the HD 142 and HD 159868 systems, further monitoring has revealed additional 
planets, and the best-fit eccentricities of the previously known planets have significantly 
decreased. Figure [8] shows the distribution of eccentricity versus semimajor axis for single 
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planets (open circles) and multiple planets (filled circles). A K-S test shows that there is only 
a 2.9% probability that the eccentricities of single- and multiple-planet systems are drawn 
from the same distribution. Previous analysis of the properties of multiple-planet systems 
have given the same result: that planets in multiple systems tend to have lower eccentricities 
than single planets (jWittenmyer et al.l 120091 ; IWright et al.l 120091 ). To investigate whether 
this difference arises from an observational bias, we asked the following question: "Is there 
a minimum threshold number of observations N for which single and multiple planets have 
the same eccentricity distribution?" We repeated the K-S tests on subsets of the known 
exoplanet data, including only those planets which have more than N observations for a 
range of N, as shown in Table [TTJ If there is a bias arising from the number of observations, 
then for larger iV & s , there would be no significant difference in the eccentricity distributions of 
single and multiple planet systems. The K-S significance levels in Table [TT] do not show any 
consistent trend with iV & s . We conclude that the difference in the eccentricity distributions 
of single and multiple planets is real and does not arise from the observational sampling. 

Both HD 142b and HD 159868c have low-eccentricity orbits near one year about stars 
similar to the Sun. This raises the question of the potential habitability (e.g. Horner & 
Jones 2010) of terrestrial moons which may orbit these giant pla nets, or Trojan companions 
of those planets, a topic explored in detail by lTinney et al.l ( 120111 ) for HD 38283b, a gas-giant 
planet in a one-year orbit. 

While it is somewhat speculative to discuss the potential habitability of as- yet undiscov- 
ered moons or planet-mass Trojan companions of planets such as HD 142b and HD 159868c, 
it is important to note that, at least for our own Solar system, the capture of objects to 
su ch orbits is now con sidered a well established part of planetary formation and migration. 
In iTinney et al.l (120111 ). we provide an extensive review of the research into satellite forma- 
tion and evolution, but it is worth reminding the reader of a few salient points. First, the 
combined mass of the largest satellites of the gas giant planets in our Solar system typically 
amounts to 2.5 x 10 -4 that of their host planet. In other words, it seems reasonable to 
expect that the most massive satellites of HD 142b could well be similar to, or somewhat 
more massive, than the Galilean satellites (since HD 142b is around 1.2 times the mass of 
Jupiter), whilst those of HD 159868c would likely be slightly less massive. As such, the 
regular satellites of those planets (assuming the same formation mechanism as the Galilean 
satellites) might well be somewhat too small to host sufficient atmosphere to allow liquid 
water on their surface. However, larger and more habitable satellites are clearly not beyond 
the bounds of possibility - particularly when one considers the possibility of the capture of 
massive irr egular satellites (such as Neptune's moon Triton) during the course of the planet's 
migration ( IJewitt &: Sheppardll2005l : IJewitt &: Haghighipourl 120071 ). 
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A more promising alternative in the search for habitable exoplanets in these systems 
could be the capture of objects as Trojans of the planets in question. Within the Solar 
system, it is known that objects can be temporarily captured as Trojans for long periods of 
time, even in the abse nce of planetary migration (e.g. Horner &: Evans 2006). However, it is 
now widely accepted (ILykawka fc Hornerll2010l ; iLykawka et al.ll2009l ; iMorbidelli et al.l 120051 ) 
that the migration of the giant planets resulted in their capturing significant populations of 
Trojans. In the case of Jupiter and Neptune, those Trojan populations were captured on 
orbits of sufficient stability that they have survived to the current day. If either of HD 142 b 
and HD 159868 c was able to capture a sufficiently large planetary embryo as a Trojan during 
their inward migration to their current location, it is highly likely that such an object could 
remain trapped as a Trojan for the lifetime of the planetary system. If such planet-mass 
Trojans exist in either system, they could well represent potentially habitable worlds. 
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Table I. AAT Radial Velocities for HD 142 



JD-2400000 Velocity (ms- 1 ) Uncertainty (ms -1 ) 



50830.95872 


29.8 


3.2 


51121.01944 


36.9 


3.6 


51385.31047 


84.7 


5.9 


51411.20252 


75.5 


6.3 


51473.08503 


33.9 


3.3 


51525.92509 


54.0 


4.0 


51526.95290 


46.3 


2.9 


51683.33138 


86.4 


3.8 


51743.27654 


85.3 


3.4 


51745.26417 


90.7 


5.5 


51767.26990 


63.9 


3.6 


51768.25417 


66.1 


3.2 


51828.06072 


40.4 


3.9 


51856.06429 


12.0 


5.5 


51856.92498 


46.3 


7.4 


51918.94072 


55.5 


3.8 


52061.29661 


84.6 


3.4 


52092.26831 


70.0 


3.3 


52093.28756 


59.4 


3.3 


52127.22295 


34.6 


4.1 


52128.15455 


42.4 


3.9 


52130.24335 


35.4 


3.6 


52151.21126 


31.6 


2.9 


52152.07857 


36.5 


3.7 


52154.15414 


28.8 


3.3 


52187.09998 


34.9 


2.9 


52188.03596 


30.0 


2.9 


52189.01990 


22.3 


3.0 


52190.00244 


29.8 


2.9 


52423.32977 


58.3 


3.1 


52425.33759 


38.8 


3.2 


52456.32088 


15.8 


3.5 


52477.24794 


25.6 


3.4 


52511.09845 


37.8 


3.4 


52654.91588 


59.7 


4.7 


52784.32722 


36.4 


3.4 


52857.23832 


24.1 


2.3 


52861.30634 


24.1 


5.1 


52946.03964 


27.6 


3.8 


53007.97796 


35.1 


3.9 


53041.95410 


72.0 


4.4 


53042.90882 


64.9 


4.9 


53215.27934 


-7.3 


2.6 


53243.28286 


6.6 


3.7 


53244.22676 


0.4 


3.8 


53246.09497 


4.5 


2.9 
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Table 1 — Continued 



JD- 2400000 


Velocity (m s 1 ) 


Uncertainty (ms : ) 
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54374.14271 


-24.4 


2.1 


54428.94253 


-3.0 


1.7 


54780.10864 


-3.4 


1.5 


55076.28334 


-58.7 


2.1 


55101.11333 


-19.8 


4.2 


55170.91268 


-7.3 


2.0 


55376.32957 


-89.9 


3.8 


55377.32745 


-86.0 


2.7 


55401.20729 


-59.7 


2.3 


55428.24962 


-16.4 


4.0 


55457.13465 


-16.3 


2.5 


55518.99948 


32.5 


2.3 


55521.00990 


10.7 


2.6 


55751.31147 


-31.3 


3.6 


55757.24366 


-34.3 


2.7 


55786.25163 


0.7 


3.2 


55788.32904 


-22.7 


3.2 


55845.10610 


14.7 


2.6 


55874.03173 


20.5 


3.0 


55874.98064 


15.4 


3.6 


55897.91075 


7.6 


2.3 
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Table 2. AAT Radial Velocities for HD 159868 



JD-2400000 Velocity (ms^ 1 ) Uncertainty (ms^ 1 ) 



52390.22780 


29.1 


1.4 


52422.14712 


10.2 


1.4 


52453.04281 


11.9 


1.5 


52456.07120 


20.2 


1.6 


52477.02060 


8.1 


1.5 


52711.26892 


-29.3 


2.9 


52747.25274 


-34.7 


1.5 


52751.26030 


-34.4 


1.4 


52786.09989 


-59.4 


1.3 


52858.95394 


-25.3 


1.3 


52942.93463 


-17.4 


1.9 


53214.09869 


8.2 


1.8 


53216.04487 


9.1 


1.4 


53242.96824 


20.6 


1.3 


53484.23874 


31.4 


1.3 


53486.16508 


35.1 


1.4 


53510.17399 


18.3 


1.4 


53521.19560 


17.8 


1.4 


53572.06973 


34.9 


1.3 


53631.89779 


49.5 


1.2 


53842.24103 


-31.5 


1.5 


53939.00914 


-16.8 


1.3 


53947.05674 


-21.3 


1.2 


54008.91766 


-9.9 


0.9 


54011.90938 


-5.1 


1.4 


54015.94916 


-9.4 


1.1 


54017.89861 


-3.0 


1.3 


54037.89603 


-9.0 


1.4 


54224.26216 


-40.6 


1.2 


54227.16135 


-28.6 


1.3 


54255.03411 


-27.9 


1.2 


54371.88740 


25.7 


1.1 


54553.20262 


18.1 


1.6 


54907.25523 


4.8 


1.7 


55101.91055 


1.3 


1.2 


55104.94555 


-10.7 


1.4 


55109.96044 


-22.7 


1.3 


55313.25667 


-52.9 


1.4 


55317.16029 


-35.3 


1.5 


55376.12410 


-16.8 


1.5 


55399.06643 


-7.1 


1.4 


55429.85887 


6.8 


1.4 


55456.87174 


4.8 


1.7 


55664.29971 


19.6 


1.3 


55692.24347 


36.7 


1.5 


55756.91450 


74.4 


1.9 



Table 2 — Continued 



JD-2400000 Velocity (ms- 1 ) Uncertainty (ms- 1 ) 
55786.07264 65.2 1.9 



Table 3. Keck Radial Velocities for HD 159868 



JD-2400000 


Velocity (m s 1 ) 


Uncertainty (ms - 1 ) 


54246.99415 


-39.3 


0.8 


54247.97421 


-32.3 


1.1 


54248.92666 


-35.4 


1.1 


54251.95192 


-24.7 


1.1 


54255.97421 


-24.2 


0.8 


54277.88461 


-17.7 


1.2 


54278.93840 


-13.8 


1.1 


54304.82235 


2.3 


0.8 


54305.83320 


6.3 


0.7 


54306.82686 


0.0 


1.0 


54307.86512 


3.6 


0.7 


54308.85011 


-0.6 


0.7 


54309.83390 


-1.1 


0.7 


54310.82821 


6.4 


0.7 


54311.82166 


9.2 


0.7 


54312.81832 


8.5 


0.7 


54313.79003 


9.7 


0.7 


54314.79189 


8.9 


0.7 


54335.72425 


20.8 


1.1 


54601.91517 


24.8 


1.1 


55024.86442 


-2.8 


0.6 


55049.83900 


-11.9 


0.8 


55052.81141 


-0.7 


0.6 


55260.15738 


-60.9 


1.0 


55369.98576 


-8.6 


0.6 


55409.84249 


-10.9 


0.6 


55462.73665 


-2.1 


0.7 


55638.13563 


-1.6 


0.8 


55665.10001 


14.3 


0.7 


55670.06654 


14.8 


0.7 


55720.04108 


48.9 


0.7 


55750.87325 


56.1 


1.3 


55825.74817 


40.1 


0.7 


55839.72241 


35.0 


1.2 
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Table 4. Stellar Parameters for HD 142 



Parameter 


Value 


Reference 


Spec. Type 


Gl IV 


Tinnev et al. (2002) 




F7 V 


Gray (2006) 


Mass (M ) 


1.15±0.10 


Tinnev et al. (2002) 




1.232tg;^ 


Takeda et al. (2007) 


Distance (pc) 


25.7±0.3 


van Leeuwen (2007) 


M v 


3.66 




Radius (R©) 


1.47±0.04 


van Belle & von Braun (2009) 




1.40±0.05 


Takeda et al. (2007) 




1.43±0.07 


Lane (1980) 


V sin i (kms -1 ) 


10.4±0.5 


Valenti & Fischer (2005) 




10.35±0.50 


Butler et al. (2006) 


logR' HK 


-4.92 


Tinnev et al. (2002) 




-4.95 


Jenkins et al. (2006) 


[Fe/H] 


0.10±0.03 


Valenti & Fischer (2005) 




0.09±0.05 


Sousa et al. (2008) 




-0.02±0.07 


Bond et al. (2006) 




0.117±0.070 


Gonzalez & Laws (2007) 




-0.02±0.06 


Ramirez et al. (2007) 


Teff (K) 


6403±65 


Sousa et al. (2008) 




6150±35 


Bond et al. (2006) 




6245±48 


Malvuto & Shvelidze (2011) 




6249 


Valenti & Fischer (2005) 




6170 


Randich et al. ( 1999) 


log g 


4.62±0.07 


Sousa et al. (2008) 






Takeda et al. (2007) 




4.19 


Valenti & Fischer (2005) 




4.2 


Randich et al. (1999) 



Table 5. Stellar Parameters for HD 159868 



Parameter Value Reference 



Spec. Type 


G5 V 


Houk (1978) 


Mass (M ) 


1 ns7+ - 032 
l-U8/_ 033 


Takeda et al. (2007) 






Valenti & Fischer (2005) 




0.919 


Sousa et al. (2008) 


Distance (pc) 


52.7±3.0 


Pcrrvman et al. (1997) 


M v 


3.63 




V sin i (kms -1 ) 


2.1 


Valenti & Fischer (2005) 


1 °S R 'hk 


-4.96 


Jenkins et al. (2006) 


[Fe/H] 


0.00 


Valenti & Fischer (2005) 


-0.08±0.01 


Sousa et al. (2008) 


Teff (K) 


5623 


Valenti & Fischer (2005) 




5558±15 


Sousa et al. (2008) 


log a 


o QQ + 0.05 

o.yy_ 04 


Takeda et al. (2007) 




3.92 


Valenti & Fischer (2005) 




3.96±0.02 


Sousa et al. (2008) 
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Table 6. Keplerian Orbital Solutions 



Planet 


Period 


To 


e 




K 


M sin i 


a 




(days) 


(JD-2400000) 




(degrees) 


(m s _1 ) 


(Mj up ) 


(AU) 


HD 142 b 


349.7±1.2 


52683±26 


0.17±0.06 


327±26 


33.2±2.5 


1.25±0.15 


1.02±0.03 


HD 142 c 


6005±477 


55954±223 


0.21±0.07 


250±20 


55.2±3.0 


5.3±0.7 


6.8±0.5 


HD 159868 b 


1178. 4±8. 8 


53435±56 


0.01±0.03 


350±171 


38.3±1.1 


2.10±0.11 


2.25±0.03 


HD 159868 c 


352.3±1.3 


53239±21 


0.15±0.05 


290±25 


20.1±1.1 


0.73±0.05 


l.OOiO.Ol 




Fig. 1. — Left panel: Periodogram of residuals for HD 142 after fitting one planet at P = 350 
days. A long-period signal is clearly present. Right panel: Periodogram of residuals for 
HD 142 after fitting two planets, at P = 350 and 6005 days. The highest remaining peak is 
at a period of 108 days, with a false-alarm probability of 5.1%. 



Table 7: The relative posterior probabilities of models Mk with k — 0, 3 Keplerian signals 
given the AAT data for HD 142. The velocity jitter was also fitted as a free parameter in 
the model. 



k P(M k \d) Jitter (ms" 



RMS 



ms 




1 

2 
3 



1.5 x 10~ 48 
1.2 x lO" 32 
1.4 x 10~ 8 
~1 



39.7 
23.7 
11.2 



39.1 
22.5 
11.2 
8.6 
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Fig. 2. — Two-planet fit for HD 142. The residuals of this fit are 11.2 ms 1 , and no further 
significant signals are present. 
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Orbital Phase 
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JD-2400000 



Fig. 3. — Left panel: Fit for HD 142b, the previously known planet with P = 350 days. Two 
cycles are shown for clarity. Right panel: Fit for HD 142c, with P = 6005 days. In both 
panels, the signal of the other planet has been removed. 



Table 8. The three-planet solution of HD 142 radial velocities. MAP estimates of the 
parameters and their 99% Bayesian credibility sets. 



Parameter 



Planet b 



Planet c 



Residual 3 



quit P [days] 351.1 [348.3, 353.8] 7900 [5500, 22200] 108.39 [107.79, 109.00] 



K [ms" 1 ] 
w [rad] 
Mq [rad] 
m v sini [M Jup ] 
a [AU] 



0.15 [0.00, 0.34] 
31.6 [26.1, 37.1] 
4.3 [3.3, 6.2] 
2.6 [1.4, 5.1] 
1.21 [0.88, 1.54] 
1.028 [0.915, 1.120] 



0.18 [0.00, 0.72] 
52.6 [40.0, 75.9] 
3.2 [1.8, 4.9] 
3.0 [0, 2tt] 
5.5 [3.7, 11.6] 
8.0 [6.0, 17.8] 



0.12 [0.00, 0.56] 
11.6 [6.1, 16.5] 
4.7 [0, 2tt] 
5.3 [0, 2tt] 
0.30 [0.15, 0.45] 
0.469 [0.418, 0.511] 



7 [ms ] 

/t [mc 1 1 



8.6 [-5.6, 74.1] 
8.8 [6.3, 11.8] 



a In the Bayesian analysis, we fit a third planet to illustrate the uncertainties in its 
parameters; we do not as yet claim a third planet in this system. 
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Fig. 4. — Periodogram of AAT residuals for HD 159868 after fitting one planet at P = 1178 
days. An additional signal is present near 355 days. 
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Fig. 5. — Left panel: Two-planet Keplerian fit for HD 159868; Circles are AAT data, and 
triangles are Keck data. The total RMS about of this fit is 5.8 ms" 1 . Right panel: Fit for the 
355-day planet only, folded to show phase coverage. The symbols have the same meaning, 
and two cycles are shown for clarity. 



Table 9: The relative posterior probabilities of models Aik with k — 0, 2 Keplerian signals 
given the AAT data for HD 159868. 



k 


P{M k \d) 


RMS 


ms 1 







4.1 x HT 23 


28 


.9 


1 


1.3 x 1(T 13 


16.0 


2 


~ 1.0 


6.6 
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Table 10. Separations for HD 142 and its Stellar Companion 



JD-2400000 


Separation (arcsec) 


Angle (degrees) 


Instrument 


Reference 


53180.7 


4.1±0.02 


184.16±0.18 


VLT-NACO 


Eeeenbereer et al. (^OO?) 


53180.8 


4.11±0.02 


184.16±0.18 


VLT-NACO 


Eeeenbereer et al. (2007) 


53316 


4.1±0.02 


184.18±0.18 


VLT-NACO 


Eeeenbereer et al. (2007) 


53712.7 


4.08±0.02 


184.13±0.29 


VLT-NACO 


Eeeenbereer et al. (2007) 


55391.43 


3.965±0.013 


184.47±0.26 


Gemini-NICI 


This work 



Table 11. Eccentricity Distributions for Planets in Single and Multiple Systems 



Filter 


K-S Probability 11 


^single 


^multiple 


All N obs 


0.029 


268 


98 


N obs > 40 


0.006 


138 


90 


N obs > 50 


0.040 


99 


88 


N obs > 60 


0.005 


64 


81 


N obs > 70 


0.088 


53 


76 


N oba > 80 


0.519 


45 


66 


N obs > 90 


0.186 


37 


62 


N obs > 100 


0.196 


30 


60 


N oba > 110 


0.082 


21 


49 


N oba > 120 


0.039 


20 


43 


N obs > 130 


0.063 


18 


43 


N oba > 140 


0.504 


12 


40 



a Probability that the two samples are drawn from the 
same distribution. 
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Fig. 6. — Dynamical stability of the HD 142 system, where planet b is fixed on its best-fit 
orbit (panel a), and a candidate planet d is present at a = 0.47AU. The color bar indicates 
the survival time in log years. Subsequent panels give the results when planet b is fixed 
on a higher eccentricity and smaller semimajor axis, in steps of la. Each square in the 
grid shows the mean lifetime of 25 independent simulations (9, for panels b and c). Panels 
(b)-(d) show the increasing instability as various mean-motion resonances (particularly the 
3:1) move with planet b. 
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Fig. 7. — Observed separations between HD 142 and its stellar compan ion. The 4 data points 
on the left hand side are from VLT-NACO (jEggenberger et al.l 120071 ). and the solitary data 
point on the right hand side is from Gemini-NICI observations. The date of observation is 
indicated next to each data point. Further details can be found in Table [TUl 
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Fig. 8. — Distribution of eccentricity versus semimajor axis for radial-velocity discovered 
planets in single systems (open circles) and multiple systems (filled circles). Planet data 
from the Exoplanet Orbit Database at exoplanets.org (2012 Jan 12). Planets in multiple 
systems are marginally less eccentric than single planets. 



